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SUMMARY 

Trajectories of water droplets about an ellipsoid of revolution with 
a fineness ratio of 5 (which often approximates the shape of an aircraft 
fuselage or missile) were computed with the aid of a differential ana- 
lyzer. Analyses of these trajectories indicate that the local concentra- 
tion of liquid water at various points about an ellipsoid in flight 
through a droplet field varies considerably and under some conditions may 
be several times the free- stream concentration. Curves of the local con- 
centration factor as a function of spatial position were obtained and are 
presented in terms of dimensionless pa-^pmeters ReQ (free-stream Reynolds 

number) and K (inertia), which contain fli^t and atmospheric conditions. 
These curves show that the local concentration factor at ’any point is 
very sensitive to change in the dimensionless' parameters ReQ and K. 

These data indicate that the expected local concentration factors should 
be considered when choosing the location of, or when determining anti- 
icing heat requirements for, water- or ice-sensitive devices that protrude 
into the stream from an aircraft fuselage or missile. Similarly, the 
concentration factor should be considered when choosing the location on 
an aircraft of instruments that measure liquid-water content or droplet- 
size distribution in the atmosphere. 


INTRODUCTION 

It is generally recognized that an aircraft moving through a cloud 
alters the concentration of cloud droplets in the immediate vicinity of 
the aircraft. For example, the peculiar distributions of ice often found 
on aircraft antenna support arms or on pitot static tubes after flight 
throu^ supercooled clouds are indications that the concentration of liq- 
uid water in the vicinity of the fuselage is considerably altered by the 
air flow about the fuselage. This effect is illustrated by the photo- 
graphs of figure 1, which show ice formed on instrument support rods 
mounted on the side and bottom of the fuselage of a B-25 aircraft during 
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flight through a supercooled cloud. As indicated by the ice formations 
of figure 1^ there is frequently a region of reduced (or zero) droplet 
concentration next to the aircraft surface, followed by a narrow region 
of greatly increased droplet concentration farther out. Beyond the narrow 
region of high concentration, the droplet concentration gradually de- 
creases toward the free-stream value with increase in distance from the 
surface of the aircraft. 

A knowledge of this spatial variation of local droplet concentration 
about an aircraft or missile during flight throu^ clouds, drizzle, or 
rain is often important when choosing the location of devices which pro- 
trude into the stream or when determining the heat required to protect 
the devices from ice. Examples of such devices are: (l) water- or ice- 

sensitive external accessories and aircraft instrument sensing elements, 
(2) Intake ducts and vents, (3) antenna masts, (4) ice detectors, and 
(5) instruments for measuring liquid-water content and droplet-size dis- 
tribution. In addition, at places where a body of revolution joins an 
airfoil (e.g., where a rocket pod is attached to a wing), the local im- 
pingement of cloud droplets on the airfoil, and therefore the heat re- 
quired for ice protection, will be altered by the effect of the body of 
revolution on the local droplet concentr-«oion. Similarly, droplet im- 
pingement on objects mounted on a wing will be affected by the variation 
in droplet concentration caused by the air flow about the wing. 

As part of an evaluation of the effect of the air-flow field on the 
droplet concentration about an aircraft or missile in fli^t through a 
droplet field, a study of the droplet concentration about a prolate el- 
lipsoid of revolution of fineness ratio 5 was undertaken. An ellipsoid 
of revolution was chosen because it is a good approximation of the shape 
of the fuselage of many aircraft and missiles. Droplet trajectories 
about the ellipsoid of fineness ratio 5 in axisymmetric , incompressible 
flow were calcxilated with the aid of a differential analyzer at the MAC A 
Lewis laboratory. Droplet trajectories were calculated as far back as 
the midpoint of the ellipsoid and as far out in the radial direction as 
1.6 times the semiminor axis. These trajectories were analyzed to deter- 
mine the relation between the droplet concentration at various points in 
space and the following variables: ellipsoid length and velocity, droplet 

size, and fli^t altitude and air temperature. The results of the analy- 
ses are summarized in this report in terms of the dimensionless parameters 
Rsq and K, which contain these variables . Although the calculations 

were made for incompressible flow, they should be applicable throughout 
the subsonic region because of the small' effect of compressibility on 
droplet trajectories (ref. l) and the hi^ flight critical Mach number of 
the ellipsoid. 
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SYMBOLS 


The following symbols are used in this report: 

A annular area perpendicular to major axis through which droplet 

flux is flowing, sq ft 

2 2 

C local concentration factor, d(rQ)/d(r ), dimensionless 

d droplet diameter, microns 

d^ed volume -median droplet diameter, microns 

F flux density of liquid water, lh/(hr)(sq ft) 

K inertia parameter, 1.704X10” d U/pL, dimensionless [the 

density of water, 1.94 slugs/cu ft, is included in the 
constant] 

K inertia parameter based on volume-median droplet diameter, 

dimens ionle s s 


L 

W 


n 

^0 

Re, 


0 


^■®0,med 


r.z 


major axis of ellipsoid, ft 
local droplet flux, number/(cm^) (sec) 
local droplet -number density, number/cc 
free-stream droplet -number density, number/cc 

free-stream Reynolds number with respect to droplet, 4.813X10~6 
dpg^U/p, dimensionless 

free-stream Reynolds number based on volume-median droplet 
diameter, dimensionless 

cylindrical coordinates, ratio to major axis, dimensionless 

ordinate of ellipsoid surface, ratio to major axis, 
dimensionless 


Tq starting ordinate at z = -» of droplet trajectory, ratio to 

major axis, dimensionless 

U free-stream velocity or flight speed, mph 

u local air velocity, ratio to free-stream velocity, dimensionless 
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u’ local air velocity, mph 

V local droplet velocity, ratio to free-stream velocity, dimensionless 

v' local droplet velocity, mph 

■w local liquid-water content, g/cu m 

w^ indicated liquid-water content, g/cu m 

Wq free-stream liquid-water content, g/cu m 

p viscosity of air, slugs/(ft) (sec) 

pg_ density of air, slugs/cu ft 

a difference hetween local air- and droplet-velocity components, 
dimensionless 

Subscripts : 

av average of quantity over area A 

r radial component 

z axial component 

0 free-stream conditions 


METHOD OF COMPUTING DROPLET TRAJECTORIES 

The droplet trajectories about the ellipsoid of revolution of fine- 
ness ratio 5 (20 percent thick) were calculated with the aid of a differ- 
ential analyzer in the same manner as those of reference 2, except that 
emphasis was placed on the behavior of the droplet trajectories in the 
space around the ellipsoid rather than on the impingement of droplets upon 
the ellipsoid surface. The trajectories were computed out to r = 0.16 
at z = 0 (fig. 2), and all nonimpinging trajectories were computed back 
as far as the midpoint of the ellipsoid (z = O) . Trajectories were cal- 
culated in the z,r plane for various values of the parameter i/K from 
0.1 to 90 and values of free-stream Reynolds number Rsq from 0 to 8192. 

In order that these dimensionless parameters have more physical sig- 
nificance in the following discussions, some typical combinations of K 
and RCq are presented in table I in terms of the length and the velocity 
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of the elllpsoidj the droplet size^ and the flight pressure altitude and 
temperature. A typical set of trajectories for i/K « 15 and Rsq = 128 

is shown in figure 3. 

The Incompress Ihle, nonviscous air-flow field used in calculating 
the droplet trajectories about the prolate ellipsoid of revolution was 
obtained from the exact solution of the Laplace equation in prolate el- 
liptic coordinates given by Lamb (ref. 3). The details of obtaining the 
velocity components in the z,r plane from Lamb's potential function are 
given in reference 2. 


RESULTS OP TRAJECTORY COMPUTATICMS 

A series of droplet trajectories about an ellipsoid of revolution 
with a fineness ratio of 5 (20 percent thick) in subsonic axisymmetric 
air flow was computed for various combinations of the dimensionless pa- 
rameters i/K and Reg- (A procedure for rapid calciilation of i/K and 

Reg from practical flight conditions is given in ref. 4.) From these 

trajectories, the variation of droplet concentration (or flux density) 
with radial distance from the ellipsoid was determined at three positions 
along the axis of the ellipsoid. 


Average Mass Flux Density of Water In Droplet Form 

The average mass flux of water in droplet form per unit area through 
an annular area of space (of width r^ - r^, (fig. 2)) perpendicular 

to the major axis of the ellipsoid is obtained from the law of conserva- 
tion of matter. Consider the liquid water in droplet form moving between 
two surfaces formed by rotating two neighboring trajectories in the r,z 
plane about the axis of the ellipsoid as shown in figure 2. Then, 

WgUAg =! Wg^^ v^^g^^ (l) 

and 

^av _ "^av^z,av _ ^ _ ^0,2 ~ ^0,1 _ ^^^0^ . , 

Fg WgU A ^2 _ ^2 A(r^) 

The subscript 0 refers to conditions at large distances ahead of 
the ellipsoid (free-stream conditions), and the subscript av refers to 
the average of a quantity over the annular area A. The annular area A 
is measured in a plane perpendicular to the major axis of the ellipsoid 



6 


NACA TN 3153 


(fig. 2). From equation (2)^ the average flxix density through an annular 
area A or a sector of the annular area A can he written as follows : 


F 


av 


0.33wqU 


A(r§) 

A(r^) 


lh/(hr) (sq ft) 


(3) 


The constant 0.33 is a conversion factor for the units used. 

2 2 

Curves of rQ as a function of r for various values of Rcq and 

i/K at constant z positions corresponding to the midpoint region 
(z =s O) , a position l/4 major axis length from the nose (z = -0,25), and 
the nose region (z = -0.5), which were obtained from the calculated 
droplet trajectories, are presented in figure 4. 


Equation (3) and figure 4 can he used, for example, to determine the 
average mass flux density of water in droplet form passing through a sec- 
tor of an annular area A such as the entrance area of a small thin- 
lipped inlet (with inlet velocity ratio eqml to l) attached to the el- 
lipsoid moving at velocity U through a droplet field with a liquid-water 
content of wq . 


Local Liquid-Water Flux Density 

The local urnss flux of water In droplet form per unit area perpen- 
dicular to the major axis at a point in the vicinity of the moving ellips- 
oid can he obtained from equation (3) by letting A(r2) approach zero. 
Then 


F = 0.33WqUC, lh/(hr)(sq ft) 


(4) 


where 

C = d(r2)/d(r2) 

The local concentration factor C at a point in space is obtained from 
the slope of the curves of r^ as a function of r^ (fig. 4) at the 
point of interest for the Re^ and K combination being considered. 

The concentration factor as a function of r for selected values 
of ROq and K at z - 0, z = -0.25, and z = -0.5 is given in 
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figure 5. Examination of figure 5 shows that the concentration factor 
approaches 1 (free- stream value) at large values of r for all z posi- 
tions. At z = -0.50^ the concentration factor decreases to values less 
than 1 for small values of r. At the z = -0,25 and z = 0 positions^ 
the concentration factor increases from values near 1 at large values of 
r to peak values at small values of r. In addition, for most values of 
Roq and K at the z = -0.25 and z = 0 positions, there is a region 

of zero concentration between the surface of the ellipsoid and the peak 
value of the concentration factor. Because of relatively low accuracy in 
obtaining values of local concentration factor greater than 4, the average 
value of 0 over an r interval of 0.001 or 0.002 units (as indicated 
by a bracket in the figure) is given instead of the peak value whenever 
the peak value is greater than 4. 

In order to determine G for various values of Roq and 1/K not 

given in figure 5, it is necessary to interpolate or extrapolate the data. 
The peculiar variation of C with Roq and i/k makes this interpolation 

or extrapolation difficult. Therefore, as an aid in interpolation and 
extrapolation, the peak value of C of figure 5 for the z = 0 and 
z = -0.25 positions is plotted as a fu-ccion of i/k for constant Roq 

in figure 6. Althou^ peak values of C greater than 4 were not pre- 
sented in figure 5 because of their relative inaccuracy, they are used in 
figure 6 in order to facilitate interpolation and extrapolation. In addi- 
tion, the r position of the peak value of C at z » 0 and z » -0.25 
is plotted as a function of 1/K for constant ROq in figure 7. With 

the use of figures 6 and 7, the peak value of C and its r position 
can be determined for the value of Reo and i/k of interest. With the 
position and value of the peak obtained by this method, the remaining 
portion of the curve of C as a function of r can bC obtained by using 
the shape of the curve corresponding to the nearest values of Re^ and 

1/K as a guide . 

The local liquid-water content in grams per cubic meter can be ob- 
tained at any point in space by dividing the droplet flux by the z- 
component of the local droplet velocity at that point: 

WqU 

w = C, g/cu m (5) 

z 

Inasmuch as the z-component of the local droplet velocity at a given 
point in the vicinity of a body of revolution is not readily determinable, 
it is usually necessary to estimate it from the local air velocity. The 
difference between the dimensionless local air- and droplet- velocity com- 
ponents as a function of radial distance r at the z =* 0 position of 
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the ellipsoid is shown in figure 8. Examination of this figure shows 
that the dimensionless local air- and droplet -velocity components about 
an ellipsoid of fineness ratio 5 are often very nearly equal at the 
z « 0 position. Examination of the air-flow field about the ellipsoid 
(ref. 2) and figure 8 shows that at z =* 0 the ratio U/v^ is usually 

nearly equal to 1 (within 6 percent) . 


Local Droplet-Number Flux Density 

By reasoning similar to that of the preceding section^ the local 
droplet-number flux density at a point in the vicinity of the ellipsoid 
moving through a droplet field can be written in the form 

N = 44.7nQUC, number/(cm^) (sec) (6) 

The constant 44.7 is a conversion factor for the units used. The 
concentration factor C is obtained from figure 5. 

If, instead of the droplet -number il^xx density, the instantaneous 
number of droplets per cubic centimeter at a point in space in the vicin- 
ity of the ellipsoid is desired, then the local flux density at the point 
must be divided by the z-component of the local droplet velocity at the 
point. That is. 


n = C, number/cc (7) 

v' 

z 

Again, since the z-coir^jonent of the local droplet velocity is usually 
unknown, it must be estimated from the local air velocity at that point, 
as discussed in the preceding section. 


Shadow Zone 

The region of zero concentration adjacent to the surface of the el- 
lipsoid, which is evident in figure 5, will be called the shadow zone. 

This region is protected from droplet penetration by the air-flow char- 
acteristics ahead of and in the vicinity of the forward positions of the 
ellipsoid. The thickness (r - rg) of the shadow zone at each z position 

of the body for various Roq and K values is given in figure 9 . Gen- 
erally, the thickness of the shadow zone increases as z approaches 0. 
The shadow- zone thickness is 0 at the nose for all values of Rbq and 

K shown in figure 9 and becomes of finite thickness at a z position 
that depends on Rcq and K. 
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Local Concentration in Clouds witL Nonuniform Droplet Size 

The data presented in figures 4 to 9 apply directly only to flights 
in clouds^ rain^ or drizzle composed of droplets that are all uniform in 
size. The droplets encountered in the atmosphere at any instant may, 
however, he of mixed sizes. Kierefore, if the free- stream droplet-size 
distribution is known or can be estimated, the data must be accordingly 
modified (or weighted) before the local flux density of liquid water at a 
point in the vicinity of the ellipsoid can be found. 

For nonuniform droplets, the local flux density can be determined 
from equations (4) or (6) by using the weighted value of the concentration 
factor C that corresponds to the droplet-size distribution present in 
the cloud. The weighted value of C can be obtained by plotting C for 
each droplet size in the distribution (based on values of 1/K and Reo 
corresponding to each droplet diam.) as a function of the cumulative 
volume (in percent) of water corresponding to each droplet size and in- 
tegrating the resultant curve. The details of this method of weighting 
are given in reference 2. 

CM 

qJj The effect of droplet-size distri'^utlon on the concentration factor 

^ is illustrated in figure 10. In this figure, the variation with r of 

the concentration factor about the ellipsoid, at the z » 0 position for 
flight through clouds with a Langmuir "C" droplet- size distribution 
(ref. 5, reproduced in table II herein), is compared with that for flight 
through clouds with uniform droplets having a diameter equal to the volume - 
median diameter of the "C" distribution in each of four cases . 

As the droplet-size distribution becomes broader, the peak value of 
the concentration factor considerably decreases and the former shadow zone 
becomes an area with a small (but increasing with r) concentration fac- 
tor, rather than one with zero concentration factor. 


DISCUSSION OF RESULTS 

The variation of local concentration factor with spatial position 
and Roq and i/K presented in figure 5 indicates that the droplet con- 
centration in the vicinity of the ellipsoid differs considerably from the 
free-stream concentration. Under some circumstances, the local droplet 
flux per unit area may vary from 0 to many times the free-stream value 
within a short distance. As shown in figure 10, this effect is greater 
for droplets of uniform size than for clouds with a distribution of drop- 
let sizes, because of the partial "averaging out" of maximum and minimum 
regions for droplets of varying sizes. In addition, the concentration 
factor is very sensitive to small changes in 1/K and Roq, as can be 



10 


MCA TN 3153 


seen from figure 5. This sensitivity to i/K and Rsq is hrought out 

more clearly in the variation of peak value of concentration factor with 
i/K and Roq presented in figure 6. Figure 6 shows that at z = -0.25 

and z = 0^ there is a narrow region of i/k for each Rcq in which the 

peak value of concentration factor is large. For values of 1/K larger 
or smaller than those in the narrow-peak region, the peak value of the 
concentration factor decreases rapidly. 

The variation of concentration of liquid-water content with z and 
r position, droplet-size distribution, and med 

dicates that care should he exercised in locating instrument sensing ele- 
ments and small inlets or vents that protrude from the surface and are 
sensitive to impinging water or ice formation. If possible, they should 
be located where a minimum concentration factor exists . When it is nec- 
essary to provide ice protection for these protuberances, consideration 
should be given to the concentration factors expected. In some cases, 
part of the device may accumulate ice at a rate which would be several 
times that expected from the free-stream liquid-water content. 

When estimating the effect of a body of revolution on droplet flux 
in its vicinity from the data calculated from the ellipsoid, the degree 
of aerodynamic and physical similarity must be -considered. The slope of 
the surface at the z position of interest, the bluntness of the nose, 
and sudden changes in slope, such as those due to windshields, are par- 
ticularly important. 


Comparison with Observed Icing Deposits 

Qualitative observations of ice formations on aircraft pitot tubes, 
antenna support arms, and other shaft-like protuberances into the air 
stream bear out the shape of the curves shown in figure 10. Inasmuch -as 
the concentration factor determines the shape of the ice formation on a 
rod of small diameter (greater than 90-percent collection efficiency) , 
these curves give directly the shape of the expected ice formation on such 
a rod located at the indicated positions . The thickness of the ice for- 
mation would depend, of course, on the velocity, liquid-water content, 
and time of exposure. The curve of figure 10(c) was calculated for values 
of I^so,med'> ^/^ed^ ^ position comparable to those existing when 

the support rods shown in figure 1 were iced. A comparison of the curve 
of figure 10(c) with the shape of the ice on the support rods shows that 
there is good qualitative agreement between the calculated concentration 
of liquid water and that indicated by an actual ice formation. 
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Effect of Local Concentration Factor on Liquid-Water-Content and 
Droplet-Size-Distrlbution Measurements 

Most atmospheric liquid-water-content and droplet-size-distrihution 
measurements have heen made from aircraft in fli^t. For mechanical rea- 
sons, the measiurements have frequently heen made rather close to the fuse- 
lage. For this reason, measurements of liquid- water content and droplet- 
w size distribution are affected hy the local concentration factor at the 

o point of measurement and may not indicate the desired free-stream values. 

Liquid- water- content measuring instruments with high collection 
efficiency . - For liquid-water measuring devices with high collection or 
sampling efficiency (over 95 percent for most conditions), which essen- 
tially measure local flux density directly, the indicated liquid-water 
content obtained by using the free-stream velocity will be high or low in 
proportion to the weighted concentration factor, as can be seen from the 
following relation: 


w^ = f/0.33U = WqC, g/cu m (s) 

If the z-component of the local air velocity u^ is used instead 

of the free-stream velocity in the calculation of liquid- water content, 
then, 

w = F/0.33 u' = w (u/u' )C, g/cu m (9) 

i z 0 z 

Equations similar in form to equation (8) are often used for liquid-water- 
content calculations, because frequently the only velocity known is the 
airspeed of the aircraft. 

Examples of Instruments of this type are the orifice-type icing-rate 
meters (refs. 6 and 7), the heated-wire liquid-water-content meter 
(ref. 8), and cloud-droplet cameras (refs. 9 and 10). The physical ar- 
rangement of these instruments is often such that the problem of instru- 
ment location on an aircraft is particularly important. The point of 
measurement of these instruments may be as close as 6 to 12 inches from 
the fuselage. Fortunately, many of the measurements of liquid-water con- 
tent that have been taken with instruments of this type have been made 
from aircraft of the order of 100 feet long at velocities of the order of 
300 mph. These conditions, along with cloud droplets of the order of 20 
microns in diameter, give Ecq of 128.7 and i/k of 165.6 at a pressure 

altitude of 15,000 feet and an air temperature of 1° F. Examination of 
figure 5 at the z = 0 position (which is comparable to the position on 
an aircraft fuselage at which the slope of the surface becomes O) indi- 
cates that, if all the droplets are 20 microns, the concentration factor 
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is 1.26 at 6 inches from the surface of the ellipsoid (r =* 0.105) and 
1.17 at 12 inches from the surface (r =» O.ll) . If a "C" distribution of 
droplet sizes is present with a volume-median droplet size of 20 microns j, 
then the variation of local concentration with r for the ellipsoid would 
he as given in figure 10(d) ^ which shows that the weighted concentration 
factor at 6 inches from the surface is 1.26 and at 12 inches is 1.14. 

Ihus, if an icing-rate meter were mounted at the z = 0 position on a 
100-foot-long ellipsoid of fineness ratio 5, it would measure a liquid- 
water content of the order of 26 percent high at 6 inches and 14 to 17 
percent high at 12 inches . In addition, inasmuch as these calculations 
are based on axisymmetric flow, any yawing motion of the ellipsoid in 
fli^t will introduce a time -dependent variation of the local concentra- 
tion factor that would be difficult to calculate. While the effect on 
the concentration factor in the preceding example was modest (of the or- 
der of 15 to 30 percent), the picture is altered considerably if, for the 
preceding example, the droplet size is 60 instead of 20 microns and the 
droplet-size distribution is uniform. Then, ReQ =* 386.1 and i/K = 18.4, 

and the concentration factor at the z = 0 position and 6 inches from the 
surface is 0, and at 12 Inches from the surface it is about 2.8. Thus, 
at the 6- inch position the indicated watei" content for these conditions 
is 0, and at the 12-inch position it is 280 percent too high. If, in- 
stead of uniform droplets, a "C" distribution (table II) of droplets with 
a volume-median size of 60 microns is present, -then the concentration 
factor at the z = 0 position at distances of 6 and 12 inches from the 
surface is 0.24 and 1.3, respectively. These examples illustrate the 
sensitivity of the concentration factor near the ellipsoid surface to 
changes in the physical environment and the importance of avoiding taking 
measurements of liquid--water content in a comparable region in the vi- 
cinity of an aircraft, as the correction of the indicated water content 
to free -stream content would be very involved. 

Rotating multicylinders . - If the liquid-water-content measuring de- 
vice has a low collection efficiency (which varies with drop size and 
fli^t condition), or depends on a difference in collection efficiency 
among its component parts, or both, as do rotating multicylinders (refs. 

11 and 12), the situation is more complicated, because the indicated water 
content will not necessarily vary directly with the weighted concentration 
factor at the point of measurement. This is true because the amount of 
droplet flux intercepted by each cylinder depends on its collection effi- 
ciency. The collection efficiency in turn depends on the droplet veloc- 
ity, which is one of the two quantities making up flux density (eq. (2)). 
In addition, reference 13 shows that a variation in droplet concentration 
over a set of rotating cylinders introduces a false relation between cyl- 
inder diameter and relative collection efficiency and thereby affects the 
droplet-size-distribution and the liquid-water-content measurement. In 
addition, there is the effect of the change in local droplet-size distri- 
bution (which is discussed in the following section) on the measurements 
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of liquid-water content "by the cylinders. Becavise rotating-multicylinder 
measurements of liquid- water content and droplet- size distribution are 
interrelated, a change in local droplet-size distribution will also alter 
the indicated liquid-water content at the point of interest. These com- 
plicating factors make the evaluation of liquid-water-content data taken 
with rotating cylinders very difficult if the cylinders are located close 
to the aircraft fuselage . 


Fortunately, most liquid-water-content data taken with rotating 
multicylinders are not too much in error from these factors, because the 
cylinders are, in practice, generally located far enough from the fuse- 
lage so that approximately free- stream conditions prevail. Therefore, 
the average concentration factor at the position of measurement gives a 
fair indication of the degree of error in measuring liquid-water content. 
In typical installations (ref. ll) , the lowest cylinder is about 16 inches 
from the aircraft surface and the top cylinder is about 34 inches from 
the surface. With the assumed 100-foot ellipsoid moving at 300 mph at 
15,000 feet with an air temperature of 1° F and a "C" distribution with a 
20-micron volume-median droplet diameter, the variation of concentration 
factor is given in figure 10(d). At 16 inches from the surface C = 1.095 
(about 10 percent high), and at 34 inchco C = 1.04 (about 4 percent 
hi^) . Thus the liquid-water- content measurements would be less than 10 
percent too high. Actually, the values of liquid-water content summarized 
in references 14 and 15 may have a somewhat smaller error due to this 
effect than indicated by this example, because many of the measurements 
were obtained by using the local airspeed rather than the free-stream 
velocity for calculating the total catch of each cylinder. However, in 
spite of the relatively good location of rotating multicylinders, ob- 
servers have reported occasional cases where the lower cylinder caught 
excessively large or small amounts of ice. This was particularly true if 
the aircraft had a tendency to yaw during the measurements and thus alter 
the position of the shadow and peak-concentration zones . 


Droplet-size-distribution measurements . - Because the local concen- 
tration factor is dependent on droplet size, the local droplet-size dis- 
tribution varies in the vicinity of a body of revolution. Suppose, for 
example, that the droplet-size distribution of a cloud is given by 

n^ = n, + n„ + n + . . . n. = droplets/cc 


where 


nj = number of droplets in the 


jth 


size group. 


Then, since each 


droplet-size group will have a different concentration factor at a point 
in the vicinity of the body of revolution, there results from equation (?) 
(assuming u/v^ nearly equal to 1, as it usually is at the z = 0 and 

z 5= -0.25 positions) 


n = niCi + ngCg + n 3 C 3 + . . . njCj 
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or 


n =ay~^ droplets/cc (lO) 

j 

Thus, both the local number of droplets per cubic centimeter and the 
droplet-size distribution will be different from the free-stream values. 
When U/v^ is not nearly eqioal to 1 (e.g., near the nose region), then 

equation (lO) has the form 

n = U ^ n^Cj/v' (11) 

J 

Equations (lO) and (ll) clearly indicate that care should be exer- 
cised in the choice of location of droplet-size-distribution measuring 
instruments, such as a cloud-droplet camera, in order to avoid regions 
in which either the concentration factor or the ratio U/v^ differs 

greatly from 1. 


CONCLUDING EEMARKB 

The data of this report are applicable in' a quantitative manner only 
to ellipsoids of revolution with a fineness ratio of 5. They also apply 
approximately in the vicinity of the nose section of a body of revolution 
that can be physically matched with the nose section of a f ineness-ratio- 
5 ellipsoid of a given length. Because many bodies of revolution of in- 
terest are different in shape from an ellipsoid, the data of this report 
are primarily useful in pointing out in a qualitative manner the type of 
variation in liquid-water concentration that might be expected in the 
vicinity of a body of revolution in flight through a droplet field. The 
analysis indicates that, for some combinations of the dimensionless pa- 
rameters Re^ and K, the value of the local concentration of liquid 

water may vary from 0 to several times the free-stream value, depending 
on the location of the point being considered. Further, the calculations 
show that the magnitude of the concentration factor at any point is very 
sensitive to change in the parameters Re^ and K. 


Lewis Plight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 26, 1954 
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TABLE II. - "G" DROPLET-SIZE DISTRIBUTION 
(ref. 5) 

Ratio of av. drop 
diam. of each group 
to volume -median 
droplet diam., 

'^/^ed 


5 

0.42 

10 

' .61 

20 

.77 

30 

1.00 

20 

1.26 

10 

1.51 


Total liquid 
water in each 
group, 
percent 


5 


1.81 
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(Ta) Front view of ice on 1- inch- thick icing- rate-meter support rods. 

Figure 1. - Continued. Ice formations on B-25 aircraft after flight at 200 mph 
through supercooled clouds. 
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(c) Ice accumulation on pressure probe. 


Figure 1. - Concluded. Ice formations on B-25 aircraft after flight at 200 mph 
through supercooled clouds. 
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Figure 3. - Droplet trajectories about ellipsoid of fineness ratio 5 for free-stream Reynolds number Rsq 
of 128 and l/K of 15. (Ordinate scale expanded 4 times that of abscissa.) 
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Figure 4. - Continued. Square of starting ordinate as function of square of droplet ordinate at constant z-posltlon 
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Figure 4. - Continued- Square of starting ordinate as function of square of droplet ordinate at constant z-posltlon. 
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Figure 4. - Continued, Square of starting ordinate as function of square of droplet ordinate at constant z-position. 
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(f) z = Oj Beo = 8X92. 

Figure 4. - Continued. Square of starting ordinate as function of square of droplet ordinate at constant z-positlon. 











Figure 4. - Continued. Square of starting ordinate as ftmctlon of square of droplet ordinate at constant z-posltion 



Figure 4. - Continued. Square of starting ordinate as function of square of droplet ordinate at constant z-posltlon 
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Figure 4. - Continued. Square of starting ordinate as function of squeire of droplet ordinate at constant z-position 
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Figure 4. - Continued. Square of starting ordinate as function of square of droplet ordinate at constant z-posltlon. 














Figure 4. - Concluded. Square of, starting ordinate as f\inction of square of droplet ordinate at constant z-posltlon 





Figure 5. - Variation of local concentration factor with radial distance r at constant axial position 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 









Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 
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Figure 5.-- Continued. Variation of local concentration factor with radial distance r at constant axifel position 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 






Figure 5. Continued. Variation of local concentration factor with radial distance r at constant axial position 
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Figure 6. - Concluded. Variation of peak value of concentration factor with l/K for constant free-stream 
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Figure 7. - Variation of peak position of concentration factor with l/K for constant free-stream Reynolds number 
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Difference betveen local ailftl air- and droplet-Toloclty cca$cmB!atB 



Badleil distance, r 


(a) Axial component. 

Figure 8. - Difference tetween local air- and droplet-velocity ocmpcnentB against radial distance 
r at z = 0 position of ellipsoid of fineness ratio 5. 
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Axial distance, z 

(b) Free-stream Reynolds number, 128. 

Figure 9. - Continued. Thickness (r - Tq) of shadow zone as function of z-positlon. 






Axial distance, z 

(c) Pree-stream Reynolds number, 512. 

Figure 9. - Continued. Thickness (r - Pg) of shadow zone as function of z-posltion. 
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Figure 9. - Continued. Thickness (r - r ) of shadow zone as function of z-posltion. 
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Flgiire 9. - Concluded. Thickness (r - r^) of shadow zone as function of z -position. 
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NACA-Langley - 7-7-54 - 1000 


feet and air tempera- 





